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Abstract 


The lateral heterogeneity and apparent anisotropy of the upper mantle are 
studied by measuring Rayleigh and Love wave phase velocities in the period 
range 100-250 sec. Spherical harmonic descriptions of the lateral heterogeneity 
are obtained for orcier and degree up Id l=m=10. Slow regions are evident at the 
East Pacific rise, northeast Africa, Tibet, Tasman sea, southwestern North Amer- 
ica and triple junctions in the Northern Atlantic and Indian oceans. Fast regions 
occur in Australia, western Pacific and the eastern Atlantic. Details which are 
not evident in previous studies include two fast regions in the central Pacific and 
the subduction zone in the Scotia Arc region. 


Inversion for azimuthal dependence showed (1) little correlation between 
the fast phase velocity directions and the plate motion vector in plate interiors, 
but (2) correlation of the fast direction with the perpendicular direction to 
trenches and ridges. Phase velocity is high when waves propagate perpendicular 
to these structures. Severe trade-offs exist between heterogeneity and azimu- 
thal dependence because of the yet unsatisfactory path coverage. Confirmation 
of the azimuthal dependence will require more complete coverage. On the other 


hand, some details of lateral heterogeneity interpretations may 


neglected azimuthal variations. 
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1. InLroduclion 

The aims of this paper arc to obtain (l) the lateral variation of phase velo- 
city of long per.od surface waves (100-250 sec) and (2) the azimuthal depen- 
dence of phase velocity. We use the spherical harmonic approach. This has the 
advantage of not requiring an a priori assumption about the nature of the 
heterogeneity. In practice, we solve For the spherical harmonic coefficients from 
1=0 to a certain upper limit, lrriax. Thus, we can not recover sharp structural 
discontinuities. Instead, we obtain a long-wavelength averaged structure. This is 
appropriate, since only long period surface waves are used. 

There are two approaches which have been used for the inversion of long 
period (>100sec) seismograms for lateral heterogeneity of the Earth. One uses 
waveform inversion and the other uses phase velocity and/or group velocity 
measurements. The former(Woodhouse, 1983; Woodr.ouse and Dziewonski, 1983) 
is probably the most direct approach and fully utilizes the information con- 
tained in the seismograms. A disadvantage of this approach is the large storage 
and computation time required. The velocity method(e.g., Nakanishi and Ander- 
son, 1982, 1983a, 1983b, 1933c) uses only the phase information of the data; the 
amplitude information is discarded. The measured velocity is the data that the 
model should satisfy. Measurement of multiplet locations(Silver and Jordan, 
1981; Masters et al., 1982) is a varient of this approach, since essentially the 
same information is obtained and used for inversion. The velocity method 
requires relatively little storage and computation time. 

There are other advantages and disadvantages associated with the two 
methods. Waveform inversion can treat higher modes easily. It is difficult to iso- 
late modes for phase velocity and group velocity measurements. For the 
retrieval of azimuthal anisotropy, the velocity measurement approach has the 
advantage, at least at present, since it is very easy to include in the inversion 
scheme. Partial derivatives with respect to anisotropy, other than transverse 
isotropy, have not been derived. These are required in the waveform inversion 
approach. 

In the present study, we attempt to obtain the azimuthal dependence of 
phase velocity. We also obtain the lateral variation with Lmax=10, which is higher 
than previous studies(e.g., lmax=6 by Nakanishi and Anderson, 1983a, b,c). After 
presentation of the method, we discuss the problems of least-squares inversion. 
We examine what the robust features are in the solutions in the singular value 
decomposition approach and then show the results for various cases. 


2. Method 


Let t be the travel time of the phase at frequency u from the source to the 
receiver and v( y? ,'i r ) be the phase velocity as a function of position ( 8,<p ) and 
azimuth 4'. We measure 'k clockwise from north. We have 


R 

= f 


ds 


r (0,p,'k) 


( 1 ) 


where the right hand side is the line integral from the source S to the receiver 
R. To first order, this is a stationary quantity and the wave path (this line 
integral) is assumed to be along the great circular path containing S and R, 

Phase velocity at {8,<p), v( Q,<p^), can be written as 

v(0,95,*) = V 0 + v h (#9 0 + v a (0,$8,*) , (2) 


where v 0 is the spherically averaged phase velocity and v^ and v a are the devia- 
tions due to heterogeneity and anisotropy. We assume 1 v h ) , |v a j«|v 0 |. Follow- 
ing Backus(1965) and Smith and Dahlen (1973), v a can be expressed as 

v a (^,93,'{ / ) = a cos 24' 4- b sin 2'k + c cos 4T + d sin 4T . (3) 


We substitute (3) in (1), expand the denominator and obtain 
^SR 


t = 


'0 

R 
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—J v h ds 
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+ — /(a cos2^ + b s^'k + c 0054 ^ + d sin4^) ds , 


where Asr is the distance between S and R. Denoting Asr - v 0 t ~ 6cp and expand- 
ing v h , a, b, c, and d in terms of spherical harmonics, we get 

R 

v 0 6(p = £ v hlm / ds 

lm S 

R R 

+ E aim / vpi&v) cos 2'k ds + £ b lm / sin 2'k ds 

lm S lm S 

R R 

+ £ c lm / WS,?) cos 4* ds + V d lm / Yf n (6,cp) sin 4'k ds , (4) 

lm S lm S 


where v h ] m , and a] m , b] m , c ]m and d Jm are spherical harmonic coefficients of v h , 
a, b, c, and d. The left hand side, v 0 8<p, is determined from the observation and 
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we solve systems of equations like (4) for v hJm , a> ni , b] m , oi m urn! dj m . Note that 
incorporation of anisotropy 2'f' and 44' increases the number of unknowns by a 
factor of five. 

For convienience in later sections, we write the system of equations (4) as 

A x = B (5) 

where 

x = ( v h 00> v h 10 &00>" , '^00>” ,C 00'” ,, ' C ^00> ,, ‘)''' > 

B = (v 0 dp )‘ , 

and the elements of A are given by line integrals in (4). 

3. Data 

Phase velocities measured for 15 earthquakes in 1980 (Table 1) from the 
records of IDA(International Deployment of Accelerometers) and GDSN(Glcbal 
Digital Seismograph Network) network, are used in this study. The single station 
method which is used for the measurement is explained in detail in Nakanishi 
and Anderson (1983c). Source parameters used for the measurement are tabu- 
lated in Table 2 (Nakanishi and Kanamori., 1983). The data set is slightly larger 
than that of Nakanishi and Anderson( 1983a, b,c). Only R2 and R3 or G2 and G3 
are used for the analysis. Rayleigh and Love waves at periods of 100 sec, 150 sec, 
200 sec and 250 sec are analysed and the number of data used in each case is 
given in Table 3. The number of data for Love waves is less than that of Rayleigh 
waves since IDA data is restricted to the vertical component. 

Path coverage 

Figs, i and 2 show the path coverage of Rayleigh and Love waves at 200 sec. 
Overall patterns are very similar at different periods. In both Figs. 1 and 2, the 
top figure shows the number of surface wave rays that go through each 20°x20° 
block. The bottom figure shows the azimuthal coverage of rays in each block. 
The azimuth of a ray can change in a block, thus a mean azimuth is calculated 
for each ray in a block and shown here. 

It is clear from the two figures that coverage of the polar regions is not 
good, but the area near polar regions is smaller for a given (20°x20°) block size. 
Path coverage per unit area is better than shown in the figures. Good path 
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coverages near the Mediterranean, Japan, South Pacific and northern ALlanLic 
ocean correspond to the locations or the antipodes of earthquakes or stations. 
Spherical harmonic coefficients are mainly controlled by the regions with good 
coverage. The results in the poorly covered region are essentially extrapolations 
and should be treated with care. 

Overall, path coverage is good, although the azimuthal coverage in some 
blocks is not. This causes a problem, i. e . trade-off, when azimuthal dependence 
of phase velocity is incorporated in the inversion. This will be discussed later. 

4. Analysis 

Solution behavior 

First, we solved the system of equations (5) by the conventional least- 
squares approach, i.e. 

x = (A t A)- 1 A t R. 

The solution x becomes very large in some cases. Phase velocity variations 
exceeded 10 %, which suggests that considerable errors are mapped to a solu- 
tion due to the near singularity of (A T A). At this point, we adopted the singular 
decomposition approach(SVD). 

In SVD, A is decomposed as 

A = UAV T , 

where U end V are orthogonal matrix and A is a diagonal matrix whose elements 
are eigenvalues of A. Then (A t A) - 1 can be written 

(A t A)" 1 = V A -2 V T . 

We eliminate some of the smaller eigenvalues in A in order to stabilize the 
solution x. The question then arises; how does the solution change as the number 
of eigenvalues in A kept in the (A T A) -1 (hereafter p) is changed. 

Fig. 3 shows the solution behaviour of the inversion of Rayleigh waves at 20Q 
sec for heterogeneity with lmax=10. The upper figure shows the solution norm 
(x) and the residual norm (*). The horizontal axis gives the number of eigen- 
values kept (p) and the increment is taken to be 2. The solution norm is defined 
as the Eulidean norm of the vector x, i.e. 
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1 1 X I ! = ~\/ E xf . 

where x T = {x 1( x 2 ,...). The residual norm is defined as the Euclidean norm of 
Ax - P. 

As p is increased, the residual norm decreases and the resolution becomes 
better while the solution norm increases, The lower figure shows the behaviour 
of some of the solutions. Each solution is normalized for plotting purposes. Some 
of the features shown in this figure are: 

(1) There are no drastic changes in the residual and solution norms as p is 
increased. 

(2) Whenever there is a small jump in residual norm and solution norm, the 
solutions change drastically (for example, see around p= 19-23, 63-65, 71-73, 
and 109-111). 

(3) Solution behaviour is especially wild from p=113 to p=12l(maximum), 
which suggests that considerable errors are being mapped into the solution 
by incorporating these smallest few eigenvalues. 

However, even the solution for p=101 showed unrealistically large peaks and 
troughs, 4-5 %, in the phase velocity variation map. Phase velocity variation of 
this order, if present, can cause severe multipathing effects(e.g., Sobel and von 
Seggern, 1978). 

Thus we face a problem of the lack of reasonable criterion for the cut-off 
level, p. This situation does not change for Love waves or for the cases where 
azimuthal dependence is incorporated. A similar situation for the geomagnetic 
problem was reported by Whaler and Gubbins(l98l). For comparison, the case 
of Love wave heterogeneity inversion with lmax=10 is showm in Fig. 4. 

Effect of the level of eigenvalue cut-off p 

With no obvious criterion for selection of p, we can ask the questions 1.) 
what are the effects of p on the phase velocity variation map and 2.) what are 
the robust features in the map. Fig. 5 shows the Rayleigh wave heterogeneity 
inversion at 200 sec with lmax=10. From top to bottom, p is 71, 81, 91 and 101. 
Contours are drawn at each 0.5 % interval. Striped regions are positive (faster 
phase velocity) region, while the patterned regions are negative (slower). Zero 
lines (average velocity) are drawn thicker than other contours. 
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Three features we can see from this figure are as follows: 

(1) The locations of peaks and troughs do not change much as p is changed. 
Note, for example, the troughs (slower regions) in northeastern Africa, 
Tibet, near Fiji and New Zealand and the East Pacific rise to California and 
the peaks (faster regions) in the western Pacific, Australia, northeastern 
Pacific and western Africa to northwestern part of Asia. These features 
seem to be quite robust. 

(2) Although the locations are robust, the total variation increases as p is 
increased. 

(3) The shape of the zero contour changes, as p is changed. 

The same conclusions can be reached from Fig. 6, which shows the results of 
Love wave heterogeneity inversions at 200 sec with lmax=:10. In this figure, p is 
47, 53, 59 and 71 from top to bottom. 


Eject of windowing 

Once the spherical harmonic coefficients are '?>■•' ained, some kind of window 
should be applied to Lhese coefficients in ord- o /oid the ringing phenomena 
associated with truncation at Imax. Effects of windowing are shown in Fig. 7. The 
case of Love waves at 200 sec with lmax=10 and p=53 is used in this figure. 

The top figure corresponds to the Hamming window, i.e. 


w(Z) = 0.54 f 0.46 cos( 


rr 


Imax 


■l ) . 


( 6 ) 


The same w(l) is applied to all coefficients within the same angular order 1. 
The other three figures correspond to the windowing by the function 


w(Z) = cos (~ 


Z 


2 imax + n 


0 ■ 


(?) 


The parameter n is 1, 2 and 4 for the second, third and fourth figure from top. 
The effect of changing n is similar to the effect of changing the level of cut-off p, 
i.e. the locations of peaks and troughs do not change much but the zero contour 
changes as different windows are applied. As higher 1 coefficients become less 
suppressed, more details emerge. For example, a positive peak in the central 
Pacific in the top three figures is broken up into two peaks in the bottom figure. 
This is done, of course, at the risk of introducing spurious ringing in the map. 


We use the window of the bottom figure for all other figures in this paper 
unless otherwise noted. 

5, Results and discussion 

In this section, we show the results of inversions with and without the 
azimuthally dependent terms. In all cases, we set the eigenvalue cut-off level p, 
such that the solution norm does not exceed 0.008. This has the effect of setting 
the maximum value of the phase velocity variations at about 2-2.5 %. As dis- 
cussed in the last section, maximum values are not robust features of a solution 
but locations of peaks and troughs are. 

lie suits for heterogeneity inversion 

We present the results for lmax=10 in this paper. Azimuthal dependence is 
not yet taken into account. The total number of parameters is (lmax+l) s = 121 
in this case. The numbers of eigenvalues kept in the solution are given in the 
first row of Table 4. Larger p means better resolution. Rayleigh waves have 
better resolution than Love waves and in both cases the best resolution is 
achieved at 200 sec. 

Fig. 8 shows the results for Love waves. The results for 100 sec, 150 sec, 200 
sec and 250 sec are given from top to bottom. There are excellent correlations 
with surface tectonics. Generally, old oceans and shields have fast phase veloci- 
ties, while ridge regions and marginal seas have slow velocities. Subduction 
zones are generally characterized by slow velocity, which is presumably the 
effect of slow velocity in the back-arc regions. 

There are gradual changes from 100 sec to 200 sec, but the map for 250 sec 
is quite different. There is, for example, a slow peak south of Africa, which does 
not exist clearly in the top three figures. Also, a fast peak occurs in the western 
Pacific and the fast peak in the northeastern Pacific disappeared. Moreover, the 
locations of the peaks have shifted. 

There is a possibility that this is due to the relatively poor fit of the model 
to the data. As shown in the first row of Table 5a, the total variance reduction 
(TVR) at 250 sec is not as good as TVR at other periods. TVR is calculated by 
TVR = (ctq — a 2 )/ ctq , where off is the variance for the spherically symmetric 
Earth and cr s is the variance for the aspherical model obtained by inversion, But 
we do not think this is likely to be the cause for this case, since comparable 
difference of TVR exists between 150 sec and 200 sec without much difference in 


locations of peaks. It is more likely to be caused by heterogeneity In deeper 
regions of the upper mantle, which is not correlated with shallower features. 

The results for Rayleigh waves are given in Fig. 9, There are quite a few simi- 
larities with the Love wave results, For example, slow velocities in the Red Sea- 
Gulf of Aden region to the triple junction in the south Indian ocean, Tibet, East 
Pacific rise to California and Tasman Sea to New Zealand. Australia and the 
eastern Atlantic are fast in both sets of maps. 

At the same time, there exist some differences. The fast regions in the 
western Pacific are shifted somewhat. The Canadian shield is fast for Love waves 
and slow for long-period Rayleigh waves. Ridges are generally slow for Love waves 
but they are not so evident on Rayleigh wave maps. Some ridge segments are, in 
fact, fast for Rayleigh waves. Subduction zones are characterized by slow velo- 
city for Love waves but some are fast for Rayleigh waves at long periods. 

In terms of path coverage, Rayleigh waves have much better coverage than 
Love waves (Figs. 1 and 2). Nakanishi and Anderson( 1983c) discussed the prob- 
lem of source depth errors on Rayleigh wave initial phase. However, TVR of Ray- 
leigh waves is larger than that of Love waves(Tebles 5a and 5b). Thus it is more 
likely that the data should be explained by differences in penetration depth or 
transverse isotropy. Differences in subduction zones, slow for Love waves and 
fast for Rayleigh waves, are probably due to the difference in penetration depth, 
since fast velocity at depth is consistent with the subduction of cold, fast litho- 
sphere. Love wave-Rayleigh wave differences at trenches and ridges can also be 
caused by transverse isotropy. This is because SH may be less than SV due, 
presumably, to ascending and descending flow in the mantle (Anderson and 
Regan, 1983), although in general SH > SV in ‘be shallow mantle (Anderson and 
Dziewonski, 1982). 

Compa ed with previous studies by Nakanishi and Anderson(l983a, b, c; 
lmax=6), the overall patterns are very similar: slow regions at the East Pacific 
Rise, northeast Africa, Tibet, Tasman Sea, southwestern No-^h America and tri- 
ple junctions in the north Atlantic and south Indian Ocean and fast regions at 
Australia, western Pacific and the Eastern Atlantic. But a few detailed features 
have emerged. Most notable is the break-up of the fast regions in the Pacific. 
This appears clearly for Love waves at 100 sec, 150 sec and 200 sec and Rayleigh 
waves at 150 sec and 200 sec. In terms of the numbers of data and the variance 
reductions (TVR), the results at 150 sec and 200 sec are most reliable (Tablas 5a 
and 5b). Thus we believe that this apparent increase in resolution is real. 
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Another feature is the slow region in the southeast Pacific, which extends into 
the Scotia Sea for short period Love waves. 

Results for azimuthal dependence 

When azimuthal dependence is included, the number of parameters 
increases considerably. P’or heterogeneity plus 2'P azimuthal inversion, there are 
three times as many parameters and for heterogeneity plus 2^ and 4^ terms, 
there are five times as many parameters. Thus Imax has to be decreased from 
10 and we mainly discuss the results for lmax=B with 2'f / dependence in this 
paper. Even for this case, we have 145 parameters. Maximum 1, imax, is taken to 
be the same for both heterogeneity and azimuthally dependent terms in (4). 
However, 1=0 is excluded from the azimuthally dependent terms. 

Figs. 10 and 11 show the results for Love and Rayleigh waves respectively 
for Lmax=G with 2^ dependence. The lines on the maps give the direction of 
fastest phase velocity at each point. The length of each line gives the difference 
of velocity between the fastest and the slowest velocity at a point, 

For both Love arid Rayleigh waves, the fastest directions do not correlate 
well with plate motion vectors. The map for Rayleigh wave at 100 sec seems to 
have the best correlation, particularly for the Pacific and Australian plates. 
Since Love waves at 100 sec period sample the shallowest parts of the Earth 
among the cases in Figs. 10 and 11, this case should have the best correlation, if 
the fast direction is the spreading direction. However, there is little correlation. 

On the other hand, the fast direction is generally nearly perpendicular to 
plate boundaries. This can be seen along the mid-Atlantic ridge, western Amer- 
ica, subduction zones in the western Pacific and the ridges in the Indian Ocean. 
This may indicate that the azimuthally dependent phase velocity is caused by 
lateral heterogeneity rather than azimuthal anisotropy. It is natural to expect 
that velocities are different depending upon whether waves propagate along or 
perpendicular to ridges and slabs. 

It is also possible that anisotropic crystals, such as olivine, are oriented or 
re crystallized by the high stresses and high temperatures at plate boundaries 
and that the fast directions correlate with the stress field. Dziewonski and 
Anderson (1983) suggest a correlation of anisotropy with stress, in any event, 
our results Indicate that anisotropy is not uniform under a given plate. Uniform 
flow is not expected when both plate drag and counterflow are taken into 
account (Hager and O'Connell, 1978; Chase, 1979). 


We need much more data before we can be fully confident of our azimuthal 
resulLs and the above discussion should be viewed in this light There is a trade- 
off between heterogeneity and anisotropy, Figs. 12a and 12b show results for 
Love waves obtained by inversion with and without the 24> azirnulhully dependent 
terms respectively. Periods are 100 sec, 150 sec, 200 see and 250 see from top 
to bottom, in both cases, lmax is taken to be 6. The cut- off level of eigenvalues 
are listed in Table 4 with other eases of inversion, lmax- *. 

Comparisons of Figs. 12a and 12b reveal that the locations of peaks and 
troughs are not so different, but definitely there exists trade-off, since the max- 
imum values in Fig. 12a are much smaller. TVR calculated for P’ig. 12a are 40.0, 
38,0, 45.9 and 36.3 % from top to bottom and those for Fig. 12b are 36.3, 30.8, 
41,6 and 33.5 %. In Table 5a and 5b, we list the variance reductions due to 
heterogeneity (HVR) for inversions with azimuthally dependent terms. HVR for 
Fig. 12a are 30.7, 29.1, 33.5 and 25.1 %■ Thus at 100 sec , for example, out of the 
total variance reduction of 4GJ' %, 30.7 % is explained by heterogeneity and the 
rest by azimuthal dependence. Care should be taken w-ith these numbers, since 
they are not quite additive. But it is clear that there exist severe trade-offs 
between heterogeneity and azimuthal dependence, 

In Table 5a and 5b, we show the results of the inversions for lmax=4 with 24> 
dependence and with 24> and 44' dependence. For Love waves at 100 sec, for 
example, TVR increases from 32.1 % to 35.2% by including 44' terms but HVR 
decreases from 23.4% to 17.8%. This feature can be seen for all cases, i.e. TVR 
increase slightly by incorporating 44' terms in the inversion but HVR decrease 
quite dramatically, in genera!, in order to retrieve high azimuthal dependence, 
good azimuthal path coverage is required. The above results suggest that the 
azimuthal coverage of our data set is not yet complete enough to recover both 
the 24' and 44' terms of the azimuthal anisotropy. 

6. Conclusions 

Using phase velocities measured for 15 earthquakes in 19B0, lateral varia- 
tion of phase veLocity and its azimuthal dependence were obtained. Because the 
matrix in the normal equation is nearly singular, singular value decomposition is 
used instead of the least-squares approach. After examining the truncation 
level of eigenvalues and the effect of windowing, the locations of peaks and 
troughs were shown to be robust but the maximum values are not. 


Inversion for heterogeneity showed two new features. The fast regions in the 
central Pacific are resolved into several smaller anomalies. A Scotia arc anomaly 
has been identified. The overall patterns, however, are similar to previous stu- 
dies with lower Imax. 

Inversion for azimuthal dependence shows little correlation between the 
fastest velocity directions and the plate motion vectors. The fastest directions 
are generally perpendicular to ridges and trenches. The phase velocity is slower 
when waves propagate along ridges and trenches. At the same time, however, 
severe trade-off between azimuthal dependence and heterogeneity was shown to 
exist due to the yet unsatlsfatory path coverage. Thus, confirmation of the 
results for azimuthal dependence should be awaited until more data are incor- 
porated into the inversion. The relation of azimuthal anisotropy to return fiov. r in 
the upper mantle will be discussed in a separate paper. 
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FIGURE CAPTIONS 


Fig. 1: The upper figure shows the number of surface wave rays in each 20 0 x2C° 
block. The lower figure shows the azimuthal coverage of rays in each 
block. This is for Rayleigh wave at 200 sec. Path coverages are similar for 
100 sec, 150 sec and 250 sec. 

Fig. 2: same as Fig. 1 except that this is for Love wave at 200 sec. 

Fig. 3: Solution behaviour of heterogeneity inversion with lmax=10. This is for 
Rayleigh waves at 200 sec. The upper figure shows the residual norm(*) 
and the solution norm(x) as a function of the number of eigenvalues kept 
(p) in the solution, they are plotted from 11 to 121 at an interval of 2. The 
lower figure shows the behaviour of the solutions. Their behaviour is wild 
for p>90 and especially so for 113^p^l21. 

Fig. 4: Same as Fig. 3 except for Love wave at 200 sec. 

Fig. 5: Effect of the level of eigenvalue cut-off, p. This is for Rayleigh wave hetero- 
geneity inversion with lmax=10 at 200 sec. From top to bottom, p is 71, 
81, 91 and 101. Locations of peaks and troughs of phase velocities do not 
change much, while the maximum values increase as p is increased. Con- 
tour interval for this and all other maps is 0.5%. 

Fig. 6: Same as Fig. 5 except for Love waves. From top to bottom, p is 47, 53, 59 
and 71. 

Fig. 7: Effect of windowing. The top figure is the result of Hamming window (equa- 
tion 6) and the rest are the results of windowing by the function in (7). 
The parameter n is 1, 2 and 4 for the second, third and fourth figure. The 
window with n=4 is u a ed for all other figures in this paper. 

Fig. 8: Results of Love wave heterogeneity inversion with lmax=10. The results at 
100 sec, 150 sec, 200 sec and 250 sec are shown from top to bottom. 

Fig. 9: Same as Fig. 8 except for Rayleigh waves. 

Fig. 10: Results of inversion with azimuthally dependent terms for Love waves. 
The 2^ terms as well as heterogeneity terms with lmax=6 are inverted, 
Plotted in this figure are the fastest direction of phase velocity at each 
point. Length of lines specify velocity difference between the fastest 
direction and the slowest (perpendicular) direction. From top to bottom, 
periods are 100 sec, 150 sec, 200 sec and 250 sec. 

Fig. 11: Same as Fig. 10 except for' Rayleigh waves. 

Fig. 12a: Phase velocity variation obtained by Love wave inversion with 2'T azimu- 
thally dependent terms. Azimuthal dependence obtained simultaneously 


is shown in Fig. 10. Ir this and the next figures (Fig. 12b), no window is 
applied. 

Fig. 12b: Results of Love wave heterogeneity inversion with lmax~8. Locations of 
peaks and troughs are similar with Fig. 12a, but their maximum values 
are much higher. This shews the existence of trade-ofT in the inversion 
between heterogeneity and azimuthal depedence. 


No. 

m 

d 

h 

m 

s 

Lat. 

Lon. 

Ms 

Region 

1 

1 

1 

16 

42 

40.0 

38.815N 

27.7B0W 

6.7 

Azores 

3 

2 

7 

10 

49 

16.0 

54.158S 

158.890E 

6.5 

Macquarie Islands 

4 

2 

23 

5 

51 

3.2 

43.530N 

146.753E 

7.0 

Kurile Islands 

7 

3 

24 

3 

59 

51.3 

52.969N 

167.670V 

6.9 

Fox Islands 

8 

6 

9 

3 

28 

18.9 

32.220N 

114.9B5W 

6.4 

Cal-Mex Border 

12 

7 

8 

23 

19 

19.8 

12.410S 

166.381E 

7.5 

Santa Cruz Island 

14 

7 

14 

16 

15 

1.7 

29.273S 

177.154V 

6.6 

Kermadec 

16 

7 

29 

3 

11 

56.3 

13.101S 

1B6.338E 

6.7 

Vanuatu Islands 

19 

10 

10 

12 

25 

23.5 

36.195N 

1.354E 

7.3 

Algeria 

21 

10 

25 

7 

0 

7.9 

L1.982S 

170.025E 

6.7 

Loyalty Island 

22 

10 

25 

11 

0 

5.1 

21.890S 

169.B53E 

7.2 

Loyalty Island 

24 

11 

B 

10 

27 

34.0 

41.117N 

124.253W 

7.2 

N. California 

26 

11 

23 

18 

34 

53.6 

40.914N 

15.366E 

6.9 

Italy 

27 

12 

17 

16 

21 

58.8 

49.479N 

129.496W 

6.B 

Vancouver Island 

28 

12 

31 

10 

32 

11.0 

46.060N 

151.453E 

6.5 

Kurile Islands 




Table 2 

Source parameters of the earthquakes used 


No. 

depth(km) 

7 

M o (ji,10 20 ) 

5 

X 

9 

1 

9.75 

17.2 

2.38 

B6.2 

3.0 

I 

CO 

o 

3 

9.75 

29.7 

1.9 

84.0 

0.0 

-70.0 

4 

43.0 

19.3 

8.31 

70.0 

89.2 

27.0 

7 

33.0 

30.1 

2.95 

60.0 

88.1 

53.3 

8 

9.75 

15.4 

0.465 

90.0 

180.0 

140.1 

12 

33.0 

51.4 

21.5 

59.0 

93.1 

170.0 

14 

43.0 

18.0 

1.38 

70.0 

82.9 

10.0 

16 

43.0 

19.1 

1.57 

54.0 

93.5 

160.0 

19 

9.75 

30.2 

4.89 

54.0 

81.8 

225.0 

21 

33.0 

38.7 

9.30 

74.0 

93.2 

143.0 

22 

33.0 

48.8 

29.2 

73.0 

88. 2 

142.0 

24 

16.0 

31.7 

10.3 

90.0 

0.0 

49.8 

26 

9.75 

44,7 

? 04 

63.0 

275.8 

-43.0 

27 

9.75 

26.2 

1.54 

90.0 

180.0 

-37.1 

28 

33.0 

27.8 

2.90 

68.0 

89.6 

2B.3 


KEY: 

T - rise time (sec) 

M 0 - seismic moment (dyne cm) 

6 - dip angle (deg) 

X - slip angle (deg) 

<p - strike, measured clockwise from north (deg) 
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Table 4 

The level of eigenvalue cut-off(p) for each case of inversion 





Love 



Rayleigh 


Imax 

Anisotropy 

100 

150 

200 

250 

100 

150 

200 

250 

10 


43 

53 

54 

53 

49 

65 

82 

74 

4 

2* 

29 

30 

28 

36 

35 

42 

52 

49 

4 

2^ + 4^ 

34 

37 

37 

43 

47 

56 

63 

02 

6 

2'T 

40 

40 

43 

46 

52 

55 

63 

63 






Table 5b 

Rayleigh waves: total variance reduction(TVR) and variance reduction 

by heterogeneity (HVR) 




100 

150 

200 

250 j 

lmax 

Anisotropy 

TVR 

HVR 

TVR 

HVR 

TVR 

HVR 

TVR 

HVR 



(%) 

(%) 

(%) 

(%) 

10 


45. B 


64.9 


66.6 


54.5 


4 

2* 

40.3 

36.3 

56.9 

38.6 

57.1 

37.B 

45.9 

18.6 

4 

2^ + 4^ 

41.3 

32.3 

60.1 

34.2 

60.1 

35.9 

47.9 

23.9 

6 

2* 

42.7 

37.7 

64.4 

52.3 

64.0 

48.4 

50.4 

35.8 


1 
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